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Quantifier 



Data members 

(required) 

Array dimensions 7, J, K 
Scalar type 

Density array (address) 
Method-choice flags 

(as needed) 
Bound_density array 
-Free_density array 
Binding_rate array 
Creation_rate array 
-IxJxK segmentation array (address) 
Segment-surface list (address) 
.Transport model (address) 



Function members 

(scalar-valued) 

Volume_in_range( [Co], [G], [segment] ) 
lntegrate( [segment | triangulated_surface] ) 
Efflux( triangulated_surface ) 

(array-valued) 

Difference( subtrahend, [segment], [divisor], [target] ) 
Multi_image( arrayjist, algorithm, [target] ) 
lmage_process( operator, [segment], [target] ) 

+ member-setting functions 

J 
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TemporalQuantifier 



Data members 
(required) 

Array dimensions /, J, K 
^Number N of times 
Times A, fz, ... , ts 
Scalar type 
Density array list 
Method-choice flags 

(as needed) 

ound density array list 
Free density array list 
inding rate array list 
Creation rate array list 
IxJxK segmentation array (static or list) 
-Segment-surface (static or list) 
Transport model (address) 



Function members 

(scalar-valued) 

Volume_in_range( [Co], [G], [segment], [times]) 
lntegrate( [segment | triangulated_surface], [times] ) 
Efflux( triangulated_surfacejist, [times] ) 

(array-valued) 

Rate_of_change( [segment], [divisor], [target], [times] ) 
Multi_image( arrayjist, algorithm, [target], [times] ) 
lmage_process( operator, [segment], [target], [times] ) 



+ member-setting functions 
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Initialize pointer on 
first triangle. 
Volume6 = 0.0 



810 



Volume += prism volume 



last triangle? 
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No 

82K ^ i 



next 
triangle 



Yes 



Return volume6/ 
6.0 
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Initialize pointer on 
first gridpoint. 
Count = 0 . 




Return count xH 
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104 



Initialize pointer on first 
volume element. 

V = 0.0 



Find 8-corner type T with respect to 
above/below threshold C 



Interpolate edge lengths relevant to T 



V += polynomial relevant to T 




Last element? 



No 
i 



next 
element 



Yes 



Return VxH 



FV 
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1020 



1030 



1050 



1100 



1140 
1141 



Initialize 'perhaps in S' 
volume element iterator. 

V = 0.0 



FV 



All corners outside segment S? 



All corners inside S? 
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Yes 



Yes 
I 



V += fractional content 



Find 8-corner type T with respect to 
condition 'above/below threshold C 



Interpolate edge lengths relevant to T 



V += polynomial relevant to T 
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Last element? 



No 
i 



next 
element 



Yes 



Return VxH 
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Input array parameters, 

range limits &, C 2 , 
optionally segment S. 



Set range 
condition 
<C 2 




Set range 
condition 
<C 2 



FV or FV 




FV or FV 



Return difference 
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Estimate dc/dt over grid, 
in array with address D 



Select a boundary B where 
c is negligible 



Fix array for free-density values, 
with address A. 



Pass D, A and B to transport 
model inversion routine. 



Return A 
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Create array e, 
initialize values to 0.0 



Triangle index i = 0 



1410 



Compute normal vector n of triangle i 



Compute mid-point w of triangle i 



Pass w to transport model 



Receive vector F from transport model 



e[N] += n. F, e[i] = n • F/(triangle area) 
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Last triangle? 



No 



index i += 1 



Yes 



Return A 
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1430 



1440 



1441 



1450 



1470 



